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ABSTRACT 



The deprojection of axisymmetric density distributions is generally 
indeterminate to within the addition of certain axisymmetric distributions 
(konus densities) that are invisible in projection. The known class of konus 
densities is expanded considerably here through the introduction of semikonus 
functions. These functions are closed with respect to multiplication in ordinary 
space, and the real part of an arbitrary polynomial of semikonus functions is 
a konus function. This property facilitates the construction of semikonus (and 
konus) functions with tailored properties, such as asymptotic forms. We also 
develop a simple technique for constructing several classes of konus distributions 
with arbitrary density profiles in the equatorial plane. 
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1. Introduction 

The deprojection of objects, such as galaxies, clusters or X-ray emission regions, is 
a common problem in astronomy. For the case of spherically symmetric objects, formal 
deprojection by means of Abel transforms is well known (Courant & Hilbert 1962). The 
general problem of deprojection of axisymmetric objects was considered by Rybicki (1986), 
who showed, by use of the Fourier slice theorem, that deprojection is indeterminate 
in principle unless the angle of inclination i = 90°. In particular, the deprojection is 
indeterminate by the addition of any function that vanishes outside a "cone of ignorance" in 
Fourier space, with its axis aligned with the symmetry axis of the object and with half-angle 
90° — i. Such distributions have been called "konus densities" by Gerhard & Binney (1995; 
hereafter GB). 

It is important to study the range of possible konus distributions to understand 
the uncertainties in the structure of observed axisymmetric objects. In their paper 
GB investigated a single konus distribution constructed from exponentials in cylindrical 
coordinates in Fourier space. In this paper we develop several methods of constructing 
konus distributions based on two ideas. 

The first idea is the related concept of a semikonus function. These are similar to 
konus functions, except they vanish in Fourier space outside a half-cone rather than a full 
cone. Given a semikonus function, it is easy to construct a konus function by adding an 
appropriate reflected function in Fourier space. This is equivalent to taking the real part 
of the function in ordinary space. The advantage of semikonus functions is that they can 
be simply combined to form new semikonus functions. Specifically, the class of semikonus 
functions is shown to be closed under multiplication in ordinary space (convolution in 
Fourier space). For example, by simply raising one known semikonus function to arbitrary 
powers, one can construct a whole set of semikonus functions. This allows the construction 
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of semikonus (and konus) functions that decrease arbitrarily rapidly at large radii. 

The second idea is a simple method for generating semikonus functions with arbitrary 
profiles in the equatorial plane of the distribution. We define a kernel for generating 
semikonus functions by the Fourier transform of a single spherical function vanishing outside 
a sphere, and then translate it up a distance k along the k z axis until it lies inside a cone. 
Weighted averages of such kernels for a range of values of k produce semikonus functions in 
terms of sine and cosine transforms of the weighting function. 

We develop the theory of semikonus functions in §2, and derive particular classes of 
solutions in §3. In §4 we discuss the implications of these results for the interpretation of 
observed axisymmetric objects. 



2. Konus and Semikonus Functions 

We adopt the notation of GB. The space density pir) of the object is related to its 
Fourier transform A(k) by 



P(r) = / G? 3 kA(k)exp(zk-r). 



The inverse of this is 

A(k) = J d 3 rp(r)exp(-ik- r). (2) 

By assumption the object is axisymmetric. Without loss of generality, the z-axis can 
be taken as the axis of symmetry, with unit vector z. In Fourier space k, consider the cone 
C(6) with axis in direction z with half-angle equal to 6. A function whose Fourier transform 
A(k) vanishes outside this cone has the property that the corresponding space density pir) 
is invisible in projection for any inclination angle i less than 90° — 9. Rybicki (1986) calls 
this the "cone of ignorance" , and GB call any such function a "konus" density. 



- 5 - 



It is very useful to sharpen the concept of a konus density. The cone C(6) consists of 
two half-cones, denoted by C + (6) and that extend in the positive and negative z 

directions, respectively. We call a distribution function whose Fourier transform vanishes 
outside one of these half-cones a semikonus function. The distinction between the semikonus 
functions defined with respect to the positive and negative half-cones is given by an 
appropriate subscript, as A + (k) or A_(k). 

The relationship between the above semikonus functions and the konus distributions 
defined by GB is simple. Any konus distribution can be represented as 

A(k) = i[A + (k) + A_(k)]. (3) 

In order that p(r) be real, the following relation must hold, 

A_(k) = A;(-k), (4) 

where * denotes complex conjugation. The corresponding result for the spatial konus 
function is 

p(r) = 9tp + (r), (5) 

where 9£/ denotes the real part of /. 

The importance of the concept of semikonus functions is that they are closed under 
convolution in Fourier space. That is, if A^(k) and A+\k) are semikonus functions with 
respect to the half-cone C + (6) } so also is their convolution 

A + (k) = A { +\k) * A { +\k) = J d 3 k'A { +\k - k')A { +\k'). (6) 

To prove this, first note that the condition that a vector k lies within C + (6) can be written 
k • z > /«|k|, where p = cos 6. Now, if two vectors ki and k 2 each lie within the cone C + (6) } 
so does their sum k = ki + k 2 , since 

k • z = ki • z + k 2 • z > p (|ki| + |k 2 |) > p\k\, (7) 
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where the triangle inequality has been used. Applying this result to Eq. (6), if k lies outside 
the cone, then at least one of the two vectors k — k' or k' must lie outside the cone, implying 
that at least one of the two functions A^(k — k') or A^(k') vanishes for any value of k'. 
Thus the integral vanishes, establishing that A + (k) is a semikonus function. Analogous 
results apply to semikonus functions with respect to the negative cone. 

The preceding result can be trivially generalized to cases where the two functions 
and may have different cone angles, 9^ and 9^ 2 \ Since a semikonus function with a 
given value of 6 is a semikonus function for any larger value of 6 as well, it follows that their 
convolution is a semikonus function with 6 = max 

It follows from the Fourier convolution theorem that if both p+\r) and p+\r) are the 
spatial representations of semikonus functions, then their ordinary product, 

Mr) = ,f>(r)-^ 2 >(r), (8) 

is a semikonus function with respect to the same half-cone. That is, semikonus spatial 
functions are closed under ordinary multiplication. This provides a simple way to construct 
new semikonus functions from other semikonus functions. In particular, applying the result 
n times to the same function p+\r) gives the new semikonus function 

pi\r)=\fi\rj\ n . (9) 



The konus distributions corresponding to Eqs. (8) and (9) are 

Kr) = 3j[A 1) (r)-A 2) (r)], (10) 



and 



^(r) = »[^(r)j • (11) 

Since semikonus spatial functions are complex, multiplication and taking the real part 
are noncommuting operations, so it is not true, for example, that /^(r) = [//^(r)]™. 
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The imaginary parts of semikonus functions are essential to the closure property under 
multiplication. 

Eq. (9) provides a simple method for constructing semikonus functions (and konus 
functions) that decay rapidly with r = |r|. For example, if p+\r) = 0(1 /r@), then 
p+^(r) = 0(l/r nf3 ). The complex phase of p^ is increased by a factor n in this process, so 
the associated konus density p^ will show a greater number of sign alterations in addition 
to the more rapid asymptotic decay. 



3. General Families of Konus and Semikonus functions 

We now shall construct several families of analytic semikonus and konus functions. 
Although they are quite simple, they provide useful examples of the kinds of indeterminacies 
that limit the deprojection of axisymmetric objects. We shall use the spherical coordinate 
r as well as the cylindrical coordinates z and R = (r 2 — z 2 ) 1 / 2 , sometimes in various 
combinations. 

Let f(k, ko) be a spherical function in Fourier space that vanishes for all k > k 0} with 
corresponding real space density 

1 r^o qiTi Jct' 

F(r,k ) = — edk——f(k,k ). (12) 
Ztt Jo kr 

We call the Fourier transform pair f(k } k ) <->■ F(r, k ) a kernel for producing a family of 
semikonus functions. If we translate a kernel along the k z axis by distance \k\ = k / sin#, 
then the non-zero region of the kernel lies entirely inside a cone of half-angle 6, and the 
density function 

p+(r,z\k,9) = F(r } k sin 9)e~ ikz (13) 

is a semikonus function that is invisible at inclinations i < 90° — 9. In general p+(r, z\k,6) is 
rapidly oscillating, but we can smooth it using the fact that any weighted sum of semikonus 
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functions is also a semikonus function. If g(k) is a weighting function, then we can build 
the smoothed semikonus function 



p+(r,z\9)= dkg(k)e- tkz F(r } k sin 9). (14) 
Jo 

By construction, this is a semikonus function for any choice of /(&, k ) and g(k). 

To facilitate evaluating the integrals in Eq. (14) we would like to have a simple form 
for the kernel function F(r, k ). The simplest kernel is a shell of radius k in Fourier space, 
f(k } k ) = 6(k — k )/k . The real space kernel is F(r } k ) = sin(A; r)/27r 2 r, and Eq. (14) 
simplifies to 

1 f°° 

p + (r,z\6) = — — / dkg(k)e~ tkz sin(kr sin 6). (15) 
2ir 2 r Jo 

Using appropriate trigonometric identities, the integral can be expressed as a sum of Fourier 
sine and cosine transformations of the weighting function g(k). If we define 

a± = r sin 6 ± z, (16) 

and the sine and cosine transforms 

f c (x) = (47T 2 )" 1 / dkg(k)cos(kx) and f s (x) = (47T 2 )" 1 / dkg(k) sin(kx) } (17) 
Jo Jo 



then 



P+ (r,z\6) = + - (18) 



Given any odd function f s (—x) = —f s (x) } the function p(r } z\6) = ?R.p + (r } z\6) is a 
konus function. The density profile of this solution in the equatorial plane (z = 0) 
is P(R) = p(R,0\6) = 2f s (R sin 6)/R. In terms of this equatorial profile, the konus 
distribution throughout all space can be written 

p(r,z\9) = -[a + P(a + ) + a_P(a)] (19) 
r 

where a± = a±/ sin 6. Unfortunately, there is no simple formula for the imaginary part of 
the semikonus density, ^sp+, in terms of the equatorial profile used to define the real part. 
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The function Qy? + (r, z\6) is also a konus function for any even function f c (—x) = f c (%) } 
but it is antisymmetric in z and therefore unsuitable by itself for representing typical 
galaxies (although it may be useful in other contexts). Knowledge of the imaginary part 
is nonetheless important here, because the closure under multiplication only holds for 
complete semikonus functions. 

Although the projection properties of the solution in Eq. (18) follow from the general 
theory presented above, it is of some interest to sketch a direct proof from the projection 
integral. This can be done by changing the variables in the projection integral to da + 
and for the functions of a + and a_ respectively. For viewing angles i > 90° — 6 the 
integration limits are asymmetric and extend from a finite value to ±oo and cancellation 
cannot occur. For viewing angles i < 90° — 6 the integration limits extend from — oo to 
+oo, and the contributions of the first and second terms in $ip+ (or ^sp+) cancel. 

We can immediately use 3?p + (r, z\6) as an example of a konus function for any odd 
function f s (x). If the defining function f s (x) asymptotically declines as r~ n , then $tp+ will 
asymptotically decline as r~ n ~ x except near the curves a± = 0. Along a curve with a± held 
constant, the function only declines as r~ x . If z = r sin <f> then holding a± = c fixed for 
r ^> 1 implies that sin <f> = =psin# ± c/r, so the region where the function declines as r -1 is 
confined to an asymptotically vanishing region near the a± = curves. A function f s (a±) 
with no extrema at finite values of a± has no asymptotic anomalies, but such functions 
cannot be used to model galaxies because they must have divergent densities as |a±| — > oo. 
The imaginary part of the solution shows the same problem. Faster asymptotic decline can 
be obtained by taking powers of this solution, but asymptotic anomalies are still prominent 
along these special curves, at least relatively. 

The nth power of the semikonus function in Eq. (18), p+(r, z\6), is also a semikonus 
function. This function declines as r~ n near the curves a± = 0, and as r~ nm away from 
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the curves if p + asymptotically declines as r~ m . Alternatively, if we chose the same 
functions but use solutions with two different cones of uncertainty 6 and 6' (6 ^ #'), then 
p + (r } z\6)p + (r } z\6') asymptotically decreases as r~ x ~ n near the two degenerate curves, 
and as r~ 2n away from the degenerate curves. The somewhat pathological asymptotic 
behavior along the a± = curves appears to be characteristic of the konus functions due to 
non-analyticity along the cone in Fourier space (the GB solution also shows this asymptotic 
property). 

A particularly simple result is found for the exponential weight function 
g(k) oc exp( — bk). In this model f c = b/ib 2 + x 2 ) and f s = x/(b 2 + x 2 ) so that 

2 sin 6 

^ = 2-2, , , . h , 2 , (20) 

r z sin v — [z + loy 



an> 



d 



^ 1+=2sm9 (b 2 + a 2 + )(b 2 + a 2 _) ^ 1+=2lSm9 (b 2 + a 2 + )(b 2 + a 2 _) ^ 

The function ?R.pi+ declines as r~ 2 at all angles away from a± = 0, while the function "^Pe+ 
declines as r~ 3 away from the lines where a± = 0. Figure 1 shows contours of the first four 
konus functions in the sequence p n (r } z\6) = ?R.pi + (r } z\6), n = 1, 2, 3, 4. Near a± = the 
solutions decline as r~ n and elsewhere as r~ 2n } and in the plane the equatorial profile is 
oc (R 2 sin 2 + b 2 )~ 2n . For comparison, Figure 2 shows the solution to Eq. (18) with the 
same equatorial profile as p 2 . This solution has only a single sign alteration instead of two. 
Figure 3 shows the angular profile of the function p\ at two different radii with r ^> b to 
illustrate the differences in the rate of convergence. Figure 4 uses the solution pE(r,z\6) to 
Eq. (18) with an exponential density profile in the equatorial plane [g(k) = 1/(1 + b 2 k 2 )]. 
The first two solutions, pEi(r } z\30°) and Pei^J^ 2;|30°) = p E1 (r } z\30°) } converge to zero 
exponentially except near a± = 0, where the convergence reverts to r -1 or r~ 2 . The last 
solution, pEi(r } z\30°)pEi(r } z\50°) which is a konus solution for 6 = 50° by the theorems of 
§2, converges exponentially at all angles. 
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Other kernels also lead to simple solutions. For example, if we define the kernel by a 
filled sphere instead of a shell, then we find the solutions 



p+{r,z\6) 
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(22) 



where the ' denotes a derivative with respect to r at fixed z (i.e. a' + = sin#). If 3?/?+ is 
to have the equatorial profile P(i?sin#), then f(x) = (x/2) uP(u)du. This class of 
solutions declines as r~ 2 along lines of constant a± instead of r~ x . If we use the weight 
function g(k) = exp(-bk) then the solution is identical to the solution p\ + found by using 
the solution for a spherical shell with the same weight function (Eq. 20) and squaring it. 
Smoothing the kernel by changing from a spherical shell to a filled sphere is (in this case) 
equivalent to the convolution of the Fourier transforms of the semikonus functions for the 
spherical shell. Polynomial kernels for /(&, k ) produce kernels F(r, k ) that are polynomials 
in k multiplied by sink r or cos k r. These lead to more general forms of Eq. (22) that 
depend on functions f(a±) and f(z) and their derivatives. We can also find analytic 
solutions for kernels including a sin ak (or cos ak) dependence, and in these solutions the 
a± "coordinates" are modified to (r ± a) sin# ± z, and the solution will include all four sign 
combinations. Spatial or parametric derivatives and integrals of semikonus functions are 
also semikonus functions. 



4. Conclusions 

Far from it being difficult to construct konus distributions, it is easy to find a variety of 
solutions with arbitrary equatorial profiles. Some of these solutions are surprisingly simple. 
Given one solution, the closure property under multiplication for semikonus functions allows 
the construction of sequences of solutions with more rapid asymptotic convergence and 
increasing numbers of zones with alternating signs. As a result, we can construct functions 
that decline arbitrarily rapidly at large radii. 
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Our construction method always produces a different rate ol asymptotic decay near 
the lines a± = r sin 6 ± z = where sin 6 is the cone opening angle. Any cut in angle 
through the density profile at a fixed large radius eventually shows a sharp pair of positive 
and negative density peaks straddling the curves a± = 0. The solution found by GB 
shows similar asymptotic properties. This pathology is probably caused by the structure 
we imposed on the integral to find analytic solutions, and it is not obviously a necessary 
property of konus functions. However, a purely numerical approach appears to be needed 
to search for konus functions with smoother asymptotic profiles, since we failed to find 
analytic solutions without these asymptotic properties. 
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Fig. 1. — Contours of the konus functions 3fy?" + (r, z|30°) for n = 1 to 4. The opening angle 
of the cone of uncertainty is 6 = 30°, and each function is normalized to have a peak value of 
unity. In projection these functions are invisible for i < 60°. The contours start at ±64% of 
the peak and there is a factor of two reduction between each contour. Positive contours are 
solid, negative contours are dashed. Note that the asymptotic decay for large radii becomes 
faster for larger values of n (oc r~ 2n ) } but there are more sign alterations. 



- 15 - 



\_ l_ l_l _) =F t I 



4.5 - 



4.0 - 



3.5- 



3.0- 




Fig. 3. — The solution to Eq. (15) with the equatorial density prohle (R 2 sin 6 2 + b 2 )~ 2 and 
# = 30°. In the equatorial plane this solution has the same density prohle as the solution 
3lpi + shown in the top right panel of Figure 1. Away from the curves a± = both solutions 
decay as r~ 4 , but near a± = this solution decays as r -1 while the p\ + solution decays as 
r~ 2 . Note also the single sign change in this solution compared to the two sign changes in 
the pl + solution. Scaling and contours are the same as in Figure 1. 
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Fig. 2. — The konus function /?i + (r, z|30°) as a function of angle (z = rsinO) for r/b = 10 
(solid) and r/b = 300 (dashed) showing the development of the peaks straddling the curves 
a± = 0. The opening angle of the cone of uncertainty is 30°. The r/b = 300 curve 
was multiplied by the ratio of the radii for the two curves (30) to emphasize the different 
asymptotic rates. 
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Fig. 4. — The top left, top right, and bottom right figures show contours of 
the exponential konus functions ?R.pEi(r } z|30°), 3tpE2('r, z\30°) = ^/^(r, z|30°), and 
^■Pei(^ } z\30°)pEi(r } z\50°). The hrst two functions are asymptotically exponential except 
near a± = where they decline as r -1 and r~ 2 respectively, while the last function is 
asymptotically exponential at all angles. The scaling and contours are the same as in Figure 
1. 



